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Abstract 

> 

On 

We discuss B — ► X s ^ decay in both constrained and unconstrained supersym- 

r — ■ metric models with explicit CP violation within the minimal flavor violation 

scheme by including tan /3-enhanced large contributions beyond the leading 
order. In this analysis, we take into account the relevant cosmological and 
collider bounds, as well as electric dipole moment constraints. In the uncon- 
Q_i' strained model, there are portions of the parameter space yielding a large 

q; CP asymmetry at leading order (LO). In these regions, we find that the CP 

phases satisfy certain sum rules, e.g., the sum of the phases of the \i parameter 
and the stop trilinear coupling centralize around tt with a width determined 
by the experimental bounds. In addition, at large values of tan/3, the sign of 
y\ ' the CP asymmetry tracks the sign of the gluino mass, and the CP asymmetry 

C3 ■ is significantly larger than the LO prediction. In the constrained minimal 

supersymmetric standard model based on minimal supergravity, we find that 
the decay rate is sensitive to the phase of the universal trilinear coupling. 
This sensitivity decreases at large values of the universal gauino mass. We 
also show that for a given set of the mass parameters, there exists a threshold 
value of the phase of the universal trilinear coupling which grows with tan (3 
and beyond which the experimental bounds are satisfied. In both supersym- 
metric scenarios, the allowed ranges of the CP phases are wide enough to have 
phenomenological consequences. 
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I. INTRODUCTION 



One of the best motivated extensions of the standard model (SM) is softly broken su- 
persymmetry (SUSY) which provides novel sources for flavor and CP violation |1[ beyond 
the Cabibbo-Kobayashi-Maskawa (CKM) picture. While flavor violation may come from 
the intergenerational entries of the soft sfermion masses and tri-linear scalar couplings, their 
phase content as well as the phases of the /i parameter and gaugino masses provide sources 
for CP violation. Flavor-changing neutral current data puts stringent bounds on flavor mix- 
ings @; therefore, such entries must be suppressed as would be the case if the same unitary 
rotation which diagonalizes the quark mass matrices also diagonalizes the squark mass ma- 
trices in flavor space. In this minimal flavor violation (MFV) scheme, which is naturally 
realized in gravity mediated supersymmetry breaking and no scale models, flavor violation 
is minimal as it is generated only by the CKM matrix. In contrast, the CP violation is not 
minimal as it can spring from both the CKM matrix which leads to flavor-changing processes 
(e.g. the CP asymmetry in B — > X^ decay) and from the flavor-blind CP phases of the 
soft SUSY-breaking masses with flavor-conserving processes (e.g. electric dipole moments 
(EDMs))i. 

The effects of the SUSY CP-violation can be manifest in several observables such as: the 
mixing of the Higgs bosons ||, EDMs [|]-[J, lepton polarization asymmetries in semileptonic 
decays ||, the formation of P-wave charmonium and bottomonium resonances [fTPf , and CP 
violation in B meson decays [|TT|^T3| and mixings [I3|. Among these observables, the most 
constraining are the EDMs [|i"5|-|17|j which bound the argument of the ft parameter to be 
< 7r/(5tan/3), leaving the other SUSY phases mostly unconstrained, in both constrained 
0-0 and unconstrained 0||0,§] supersymmetry. The constraint on the phase of the \i 
parameter is lifted if the sfermions of the first two generations weigh ^> TeV as in the 
effective SUSY scenario [113] though the EDMs are regenerated at the two-loop level, and 



can still compete with the experimental bounds |L9| in certain portions of the parameter 
space. 

The rare radiative inclusive B meson decay, B — > X s 7 provides a powerful test of the 
standard model (SM) and "new physics" such as supersymmetry. The drive to reduce 
theoretical uncertainties has led to the computation of next-to-leading order (NLO) terms 

m 



in the SM 



and two-doublet models 



and a partial computation to the same order 



in supersymmetry (SUSY) 
ALEPH 



2,p3|. The measurements of the branching ratio at CLEO p4 



25 1 and BELLE [^(J give the combined result[| 



BR (B ->■ X s -i) = (3.11 ± 0.42 ± 0.21) x 10 -4 (1) 

whose agreement with the next-to-leading order (NLO) standard model (SM) prediction 



2Tj 



BR(B^X sl ) 



SM 



(3.29 ±0.33) x 10" 



(2) 



1 The combined result is based in part on an updated value from CLEO of 3.03 x 10 4 . We thank 
G. Ganis and E. Thorndike for bringing this value to our attention. 



is manifest though the inclusion of the nonperturbative effects can modify the result slightly 
27 . That the experimental result (jl]) and the SM prediction (|]) are in good agreement 



shows that the "new physics" should lie well above the electroweak scale unless certain 
cancellations occur. In addition to the branching ratio, the recent measurement of the CP 



asymmetry has been updated to |28| 



A 



BR (b -f si) - BR ( b 



CP 



s 7 



BR(6-^s 7 )+BR(6 



.s'7 



-0.079 ± 0.108 ± 0.022)(1.0 ± 0.03) 



implying a 95% CL range of 



0.30 < A C p(5^X s 7) < 0.14 



(3) 



(4) 
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In the SM, this value is calculated to be rather small: Aqp ~ 1% 

The LEP era has ended with a clear preference to large values of tan/? ||29|| , for which it is 
known that there are tan /3-enhanced SUSY threshold corrections []30| which (i) significantly 
modify |3T|] the leading order (LO) Wilson coefficients [32|, and can (ii) dominate the NLO 
contributions pT|22| . In the absence of CP violation, weak-scale SUSY satisfies the B — > 
X s i constraints at large tan /3 more easily if the /z parameter and the trilinear soft masses 
are of opposite sign [|33"| , |3"4"f . 

In this work we will analyze B — > Xgj decay in the minimal supersymmetric extension 
of the SM in connection with the CP violating phases. Above the electroweak breaking 
scale, the SUSY Lagrangian possesses two global symmetries: a continious i?-symmetry, 
U(1)r, and a Peccei-Quinn symmetry, U(1)pq. The U(1)r symmetry is broken by the 
/i parameter, the trilinear couplings Af and the gaugino masses. The U(1)pq symmetry, 
however, is sensitive to /i parameter and Higgs bilinear mass parameter B only. Treating 
the soft masses (and \i parameter) as spurions it is easy to see that theory possesses a full 
£/(l)ij x U(1)pq invariance above the electroweak scale thereby allowing for the elimination 
of two dynamical phases. Though the electroweak breaking leaves only one symmetry to use, 
the phase of the Higgs bilinear soft mass B can be always eliminated by rephasing the Higgs 
vacuum expectation values, and therefore, one more phase can be still eliminated. Depending 
on the specific structure of the soft breaking sector this invariance allows for eliminating one 
or more phases. In the constrained minimal supersymmetric standard model (CMSSM) the 
gaugino masses are universal, and therefore, all of them can be chosen real leaving the /i 
parameter and the trilinear couplings as the only sources of the supersymmetric CP violation. 
In the unconstrained model, however, one is left with more phases to generate CP violation. 
In what follows we will assume a universal phase for the masses of the SU(2) and U(l)y 
gauginos, and measure the rest of the phases (the phase of the fi parameter, the trilinear 
couplings, and the gluino mass) with respect to them. In addition, the sign conventions (e.g. 
the sign of the fi parameter relative to the trilinear couplings) for soft masses can be fixed 
from the sparticle mass matrices listed in the appendices. 

Concerning the calculational precison, we consider the B — > Xgj decay using the LO 



Wilson coefficients [32] and by incorporating beyond leading order (BLO) tan/?-enhanced 
SUSY threshold corrections |JT| in the presence of soft SUSY phases |]T| (which we hereafter 
call the BLO scheme to differentiate it from one which includes NLO QCD corrections). It 
is important to note that the LO direct CP violation is small in both constrained |35| and 



unconstrained 

0-UH- 



3~6fl SUSY models once the EDM constraints JL5-17J are taken into account 



Sec. Ill is devoted to the study of B — > X s j decay in an unconstrained supersymmetric 
model, i.e., in a model in which soft masses are not subject to the constraints from supergrav- 
ity or stringy boundary conditions at ultra high energies. After determining an appropriate 
region of the parameter space that (i) suppresses the EDMs and (ii) enhances the SUSY 
threshold corrections, the allowed ranges of the SUSY phases as well as the resulting CP 
asymmetry will be numerically estimated. In Sec. IV, we will perform a detailed analysis of 
b — » Xsj decay in the constrained MSSM with explicit CP violation by taking into account 
the tan /3-enhanced SUSY threshold corrections at the weak scale. Since the flavor mixings 
in the CMSSM are determined by the CKM matrix, the dominant contribution to the decay 
comes from the chargino and charged Higgs exchanges with flavor-blind phases playing the 



main role | 35fl . Our conclusion are given in section V. 



II. GENERAL FORMALISM WITH CP- VIOLATING PHASES 



In this section we study the inclusive radiative 5-meson decay in a general low energy 
SUSY model with nonvanishing soft phases. To a very good approximation, B — > X s ^ is 
well approximated [BO] by the partonic decay b — > 57 whose analysis can be performed via 
the effective hamiltonian 
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(5) 



Kl s K tb} K is the CKM matrix, and the operators Oi(Q) are defined in |20|| . Also 



where \ t 

listed here are the values of the Wilson coefficients at the weak scale where the electric and 
chromoelectric dipole coefficients C-j^{Qw) are decomposed in terms of the W ± , H ± and x ± 
contributions as is appropriate in the MFV scheme. 

The negative Higgs searches at LEP prefer those regions of the SUSY parameter space 



with tan/3 > 3.5 [29|, and therefore, it is necessary to improve the leading order analysis 
by incorporating those SUSY corrections which grow with tan/3. Indeed, such non- 



logarithmic threshold corrections significantly modify tree level Higgs and chargino couplings 
30|j3l| . The effective lagrangian describing the interactions of quarks with W ± , H ± and the 
charged Goldstone boson (G 1 * 1 ) is given by 



C = 92 (Qw) [K tq t L 7 M q L W+ + h. c 



+ 



92 {QW ^ ^ {QW) [K ts ^{cot(3(l + e ts tan/3) H+ + G+} 



V2M W 

92(Qw)m;(Qw) 

V2M w (l + e b * b tan/3) 



s L + h. c. 



K tb t L {tan(3 H + - (1 + e tb tan/3) G + }b R + h. 



(6) 



where the dimensionless complex coefficients e bb , e ts and e tb represent the SUSY threshold 
corrections at the associated vertices 
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where the sfermions of first two generations are assigned an average mass of Qu- The mixing 
matrices of charginos Cl,r, neutralinos Co and squarks Grr are all defined in Appendix A, 
and the loop function H. is given in Appendix B. In computing fl7|), only the gluino and 
Higgsino exchanges are considered as they dominate over those of the electroweak gauginos 
since «2,i = dii/i.^ 71 ) ^ a s = g 2 /^ 71 ): «t = \ht\ 2 /{^)- A remarkable feature of these vertex 
corrections is that they assume nonvanishing values when all SUSY masses are of equal size 



tbb -> — e 
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8tt 
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where the two independent phases are given by Q\ = Arg [/j] + Arg[m~] and 9 2 = Arg \p] + 
Arg [At]. Numerically, |e&&| ~ |e«,| ~ |e ts | ~ 10~ 2 , so that the radiative corrections in @ will 
be 0(1) when tan/5 ~ 10 2 . 

Using the effective lagrangian (§) it is straightforward to compute the W ± and H^ 
contributions to the electric and chromoelectric dipole coefficients in (EJ) 
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(9) 
(10) 



where the loop functions are given in Appendix B. These Wilson coefficients now possess 
a CP-violating character via the tan /3-enhanced SUSY threshold corrections so that they 
can play a role in the CP asymmetry. Indeed, one may treat the BLO piece in Cf r 8 (Qw) as 
the radiative corection to the CKM factor X t (so that the W boson contribution reduces to 
its LO form) but this does not eliminate the BLO corrections (and the BLO sources of CP 
violation) as the same quantity appears this time in C^ 8 (Qw) and the chargino contribution. 
Even at the LO level, the complex part of the chargino contribution grows linearly 
with tan/3 at large tan/3 [|52]j . For instance, the heavy top squark and charginos give the 
contribution 
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where the vertex factors are defined by 
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Tl|) are defined at the scale Q s (Q s ^> Qw) 



with 7 L ^ = (r L J J r^. The dipole coefficients 

which may be identified with the masses of the heavy stop or gluino |3l| , |22| . 

In the absence of the SUSY CP phases, the theoretical estimate of b — > s'j is subject 
to experimental limits on the branching ratio ([!]) as well as the bounds on the sparticle 
masses from the direct searches. When the CP phases are switched on, however, induction 
of the EDMs is unavoidable and imposes severe constraints on the parameter space. In 
super symmetric models with explicit CP violation, the EDM of a fundamental fermion / 
(first or second generation leptons or quarks) can receive contributions from both one- and 
two-loop quantum effects: 
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where rrif is the mass of the fermion. The arguments of the one-loop contribution corre- 
sponds to the gluino, chargino and neutralino exchanges, respectively. It is not surprising 
that this one-loop term [R] behaves roughly as 
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that is, the tan /3-enhanced CP-violating contributions to the Wilson coefficients are directly 
suppressed by the one-loop EDMs unless either (i) one chooses Qu large enough ||, or (ii) 
invokes a cancellation mechanism among different SUSY contributions. In fact, studies of 
both unconstrained @J^,0,§ and constrained ||-f| supersymmetry show that the one-loop 
EDMs are sufficiently suppressed if \6^\ < 7r/(5tan/3) with no constraints on other soft 
phases. In the next section, we will follow the first option, that is, we suppress the one-loop 
EDMs by taking a large enough Q 12 , say, Q12 ^ 4 TeV ||. In Sec. Ill we will discuss b — > sj 
in the constrained MSSM in which the one-loop EDMs already agree with the bounds when 
6^ is close to a CP-conserving point. 

The suppression of the one-loop EDMs, however, is not necessarily sufficient, as there 
exist two-loop contributions [|l9j which are generated by third generation sfermions. Par- 



ticularly for down-type fermions, the two-loop EDMs grow linearly with tan/? 
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where / = e or d, and use has been made of the relation vn^j cos (3 
the two-loop functions f^~ are defined by 
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where the loop function F [x] is given in Appendix B. This two-loop contribution is pro- 
portional to tan/3/M|, and can be sizable at large tan/3 when the charged Higgs boson 
is relatively light. Although one can partially cancel (|U|) by choosing the sbottom sector 
parameters appropriately, e.g., \A b \ ~ |/x|, At ~ Q& h and &j~ Or at a specific value of tan/?, 
in general, b — > sj must be analyzed in conjunction with two-loop EDMs in determining 
the allowed portions of the SUSY parameter space. 



III. B -> X s j IN AN EFFECTIVE SUSY MODEL 



In this section we will discuss the b — ► 57 decay in the framework of an effective SUSY [IS 
scenario. We choose the first two generations of sfermions to be heavy enough to suppress 
their contributions to EDMs. Clearly, for Q 12 3> Q s the contributions of the first and second 
generations to b — > sj are also suppressed |56] . 



Given the precise fit to the electroweak observables, either both stops must weigh 
(9(1) TeV or only a predominantly right-handed stop can be allowed to weigh as light 
as ~ Qw- in other words, the stop mixing angle should be small enough to have 
t 2 = —sm9re lJ ttL + cos9^t R « £#, and t\ = sin 6^e~ %r ti R + cos 6^ ~ t^. Therefore, a 
light sparticle spectrum with hierarchy Mh, M^ 2 , M x ± ~ Qw, M^ ~ m~ ~ Q s enhances the 
SUSY contribution to b — »■ S7 without spoiling the electroweak precision data p2| , |3l| . Given 
that the SM result (0) is in good agreement with the experimental result ([T]) then it is clear 
that the contributions of the charged Higgs (which is of the same sign as the SM result) and 
chargino-stop loop must largely cancel so as to agree with experiment. 

In the limit of degenerate soft masses, the tan /3-enhanced vertex corrections assume 
nonvanishing values in (|8J). In this limiting case, the PU-contribution C^ 8 (Qw) reduces to 
its LO form as the radiative corrections cancel. However, in the very same limit, the charged 
Higgs contribution C^ 8 (Qw) maintains an explicit dependence on the SUSY phases via e&6 
and e ts . Therefore, unlike the LO C^ 8 {Qw), the BLO charged Higgs contribution obtains a 
CP-violating character, and together with the chargino contribution (which contributes to 
CP violation at the LO level) they form the two key contributions to the CP asymmetry 
in the decay. The vertex correction factors ebb, ttb and e ts are of (9(10~ 2 ) in the decoupling 
limit, and they induce large corrections at large tan/? [ |31[ . 



Next one observes that in the limiting case of degenerate soft masses (H), the radiative cor- 
rections to C^ 8 (Qw) vanish identically in accord with the decoupling theorem. However, in 
the very same limit, the charged Higgs contribution C^ 8 {Qw) still has an explicit dependence 
on the SUSY phases via e&& and e ts . Therefore, unlike the LO case, the charged Higgs-induced 
BLO dipole coefficients acquire a CP-violating potential. The tan /^-enhanced vertex cor- 
rections €bb, €tb and e ts (at least their dominant pieces proportional to a s ) assume a value 
of ~ 10 -2 in the decoupling limit. Actually, they remain close to this value in most of the 
SUSY parameter space. 

Significant stop mass splitting implies that large logarithms appear when the Wilson 
coefficients (|l|) are evolved from Q ~ Q s to Q ~ Qw The chargino contribution at the 
electroweak scale Qw is therefore obtained after resummation of such logarithms pTI]] 
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v a s (Qw)/ 
where the CP-violating parts proportional to F L ' R are defined via tilded vertex factors 

7& = 7gi+ (f?)*fg. (19) 



with f Lii ; are defined by 
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Obviously, the evolution from Q s to Qw depends on the light colored particle spectrum: 
when all colored sparticles are heavy /3q = 7, when the light stop is significantly lighter than 
Q s Po = 41/6, and when the light stop, sbottoms and gluino are all much lighter than Q s 
A) = 9/2. 

After computing the IV-boson (|), charged Higgs (|10|), and the chargino (0)111) contri- 
butions to the Wilson coefficients at the weak scale, one can use the standard QCD RGEs for 
obtaining the Wilson coefficients at the hadronic mass scale Qb ~ W&. Then the branching 
fraction [DJ and the CP asymmetry [27] can be estimated directly. In the calculations below 
we will take S = 0.9, where S is a parameter determined by the condition that the photon 
energy is above a given threshold E 1 > (1/2) (1 — 5) raj,. 

The Standard Model CP asymmetry in the b — > s'y decay is ~ 1%, and therefore, this 
quantity can be quite sensitive to new physics contributions |27|| . In an effective supersym- 
metric model with LO Wilson coefficients, the CP asymmetry can be as large as ~ 8% |36| 



when the charged Higgs, charginos and the lighter stop all have masses of order the weak 
scale. In what follows, we will determine (i) the allowed rages of the SUSY CP phases, (ii) 
the size of the CP asymmetry, and (Hi) the correlation between the asymmetry and the 
branching ratio. The numerical predictions made via the exact expressions of the threshold 
corrections in (|7j) and via the limiting expressions (§) are similar to each other. To illustrate 
the effects of the threshold corrections on the branching ratio and the CP asymmetry we fix 
the values of the parameters as M~ t2 = Mh = 250 GeV, \/j,\ = \Ab\ = 150 GeV, 9j= 7r/20, 
@A b = @A t , rtio = (Mj) 2 = (1.2 TeV) 2 , and m 2 ^ = (1 TeV) 2 where the last two masses are 
needed to fix the sbottom sector for evaluating the two-loop EDMs fllSD. When evaluating 
the threshold correction (0) we consider their limiting forms (H). 

Given that the W-boson contribution alone is already consistent with the experiment, 
and that C^iQw) an d C^ 8 (Qw) have the same sign, it is clear that some cancellation is 
needed between the H ± and x contributions. In the large tan/? regime, the difference 
between the total SUSY prediction C 7 {Qw) and the experimental result C 7 xp ~ C^ (Qw) 
behaves roughly as 
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where the symbols [■ • ■] and (■ ■ •) stand for positive numerical coefficients with and without 
the loop suppression, respectively. Shown here is only a rough estimate of the actual co- 
efficients (PHl8|) after neglecting several subleading terms. For small or moderate values of 
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FIG. 1. The allowed region in the 9a~0^, plane with LO for 10 < tan/3 < 60. 
region when BLO Wilson coefficients are included is similar. 



The allowed 



tan/3, a cancellation between the chargino and charged Higgs contributions is needed, and 
this happens when 



Ou + o At ~> 



IX 



(22) 



At higher values of tan/3, threshold corrections become important, and a suppression of such 
terms occurs when 



6^ + 6-^0 



(23) 



thus imposing a constraint on the gluino phase. These rough estimates on the allowed ranges 
of the SUSY phases are confirmed in Figure [l], which shows the allowed region in 6^-6^ plane 
for the LO approximation. The allowed region when BLO Wilson coefficients are included 
is similar. Here the BLO region is valid so long as the gluino mass is positive, 9~ = 0. When 
the sign of the gluino mass is inverted, 9~ = n, the allowed region gets reflected with respect 
to the 9n = 6a line. These direct estimates are in agreement with the rough estimates in 



(^) and (g3|) above. 

In Figure |2|, we show the variation of the CP asymmetry with tan/3 for LO (upper 
window) and BLO (lower window) Wilson coefficients. The individual points correspond to 
a scan over values of 6^ and 9a taken from Figure [l|. Since the allowed values of the phases 
in the LO and BLO approximations are not the same, the surviving points in the scans will 
also differ as seen in comparing the two panels of Figure |2|. At LO precision, Aqp ranges 
from —8% to 8% uniformly, that is, for a given value of tan/3, there is no strong preference 
to positive or negative values. When BLO Wilson coefficients are included however, the 
dependence of the CP asymmetry is modified. We find that there is an observable preference 
to to positive values of A^p distorting the uniformity of the LO behavior. This particular 



9 



behaviour of the CP asymmetry results from our choice of a positive gluino mass, 6~ = 0. 
For 8~ = ir, the graph is approximately inverted (relative to that shown in the lower panel of 
Figure 0), i.e., with a preference to negative values. Therefore, the sign of the CP asymmetry 
tracks the phase of the gluino mass for most of the parameter space at large values of tan/3. 
In general, there are two main reasons that the CP asymmetry is enhanced: (i) it is 
maximized when the branching ratio is minimized, and (ii) it can be maximized due to 
specific relations among the Wilson coefficients depending on the underlying model. Indeed, 



it is known that the asymmetry in the decay behaves as ~ 10% x |C 8 (mb)| / \C 7 (m b )\ ||i"3| , |37 
and reaches the 10% level when the chromoelectric coefficient has a size similar to the electric 
coefficient, as was first pointed out in the framework of two-doublet models |37|. Depicted 
in Figure |] is the variation of the CP asymmetry with the branching ratio of the decay, 
where the enhancement of the asymmetry with decreasing branching ratio is manifest. In 
this figure, we assume m~ > 0. For most of the points plotted, there is a correlation between 
the signs of the CP asymmetry and the gluino mass. 

In this section we have analyzed the rare radiative B meson decay B — > X s 7 in the 
framework of an effective supersymmetric model ||18[ in which the sfermions in first two 
generations are heavy enough to suppress the one-loop EDMs J|,|],[7|||, and the two-loop 
EDMs |TJ| are suppressed by appropriately tuning the stop and sbottom sector parameters 
in addition to choosing small stop mixing angles (which are also required by the electroweak 
precision data). 

In this framework we find that, in regions of the SUSY parameter space where there is 
a large CP asymmetry at LO, the BLO effects induce asymmetries which are twice as large 
as the LO asymmetry at sufficiently large values of tan/3. Moreover, certain combinations 
of the SUSY soft phases are constrained to reside close to CP-conserving points in order 
to agree with experimental bounds. The size of the allowed ranges of the SUSY phases, 
which is mainly dictated by the experimental bounds on the branching fraction and the 
EDMs, is wide enough to have an enhanced production of the P-wave bottomonia in lepton 



colliders fLO] , to have observable CP violation in the Higgs system || in next generation 
of colliders, and to have observable leptonic polarization asymmetries in 5-meson decays 
||. Neither of these phenomena have been observed yet but they are all interrelated. First 
of all, the P-wave bottomonium production rate is determined by Im[C 7 ] [|K|] which has a 
richer structure in the BLO case. Secondly, in the Higgs system, where the mixing of the 
opposite CP Higgs bosons depend on sin(6 l 2 ) sm(8~), there exists sizable CP violation in the 
parameter space of Fig. 1 when 9 2 differs from ir. Moreoever, altough 9? is relatively small 
(as required by the precision data and EDM constraints) the CP-violating Higgs mixings 
are still important with relatively light Higgs bosons [^| as allowed by b — > 57. Finally, 
a simultaneous measurement of the lepton polarization asymmetry and CP asymmetry in 
B — > K*£ + £~ decay provides a good consistencey check of the 'new physics' contributions 
where both depend on the imaginary parts of the dipole coefficients C?^ (and other relevant 
coefficients), and they are necessarily enhanced in the BLO case 0. 

IV. B -» X57 IN THE CMSSM 

In this section, we will analyze the implications of the SUSY CP violating phases on 
b — ► s'j in the framework of the minimal supergravity model, or equivalently, the constrained 
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window) Wilson coefficients for 10 < tan/3 < 60, and < (6u,0a) < ^ '■ Individual points corre- 
spond to values of the two CP violating phases taken from Figure H. For tan/3 < 35, for which 
the radiative corrections are small, both asymmetries behave similarly, and fall around 2% when 
tan/3 ~ 10. The BLO CP asymmetry shown here is computed for a positive gluino mass, and to a 



good approximation, A(jp -^ — ^CP 



as m~ 



-m~ at large tan/3. 
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for m~ > 0. For most of the points shown, the sign of the CP asymmetry is correlated to the sign 



MSSM or CMSSM |[38| . In the CMSSM, universal gaugino masses mi/ 2 , scalar masses m 



(including those of the Higgs multiplets) and trilinear supersymmetry breaking parameters 
Aq are input at the supersymmetric grand unification scale. In this framework, the Higgs 
mixing parameter \x can be derived (up to its phase which does not run) from the other 
MSSM parameters by imposing the electroweak vacuum conditions for any given value of 
tan/3. In the CMSSM, there are only two physical phases to consider. Thus, given the set 
of input parameters determined by {mi/2, m Q , \Ao\, tan/3, 6^ and 9a}, the entire spectrum of 
sparticles can be derived. 

In CMSSM, the low energy mass spectrum is completely controlled by the GUT-scale 
parameters above, and one cannot decouple the first two generations of sfermions (though 
they are nearly degenerate to an excellent approximation); therefore, the Wilson coefficients 
at the SUSY-breaking scale (|PT| ) must be updated by taking Qi 2 ~ Q s . Consequently, the 
contribution of squarks in the first two generations give 



CUQs) = \cUQs) in Eq. (11) 
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(24) 



where Tl.r, is obtained by replacing the stop mixing matrix C~ with the identity matrix in 
(|T2|). If the CMSSM spectrum admits the lighter stop to be as light as ~ Qw (which can 
occur when large values of A are assumed), then the analysis of b — ► 57 proceeds as in the 
last section, in particular, the chargino contribution at the weak scale is given by ( |TTD and 
([[ID where C* 8 (Q S ) is now given by (]2T 
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The remaining colored sparticles typically have masses around the SUSY breaking scale, 
Q s . For such a hierarchy of the masses, the analysis of the Wilson coefficients differs from 
the previous cases, in particular, the chargino contribution at the SUSY breaking scale is 
now given by 
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which can be rescaled to the electroweak scale via QCD running 



01 

M 2 , 



(25) 



C 7 (Qw) 

^8 (Qw) 




cUQs) + 



^8 \Qs 





cUQs) 



(26) 

(27) 



where /3 = 7, as the colored spectrum below Q s is just that of the SM. The formulae (|2~4l , |25"D 
give possible changes in the effective theory at the weak scale in analyzing b — > 57, which 
will be numerically analyzed by taking into account cosmological and collider constraints. 

Due to the restricted parameter set, and the strong correlation among the masses of the 
sparticles, there are a number of experimental constraints which must be considered. The 
most important of these are provided by LEP searches for sparticles and Higgs bosons [39| . 
the latter constraining the sparticle spectrum indirectly via radiative corrections. The kine- 
matic reach for charginos was m x ± = 104 GeV, and the LEP limit is generally close to this 
value, within the CMSSM framework. The LEP searches for sleptons impose mg > 97 GeV, 
m^ > 94 GeV and m f 



> 80 GeV for m x < 80 GeV. Other important sparticle constraints 



are those on stop squarks t: 



nil 



> 94 GeV for m x < 80 GeV from LEP, and m~ t > 115 GeV 



for m x < 50 GeV from the Fermilab Tevatron collider 
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The lower limit on the mass of a Standard Model Higgs boson imposed by the combined 
LEP experiments is 113.5 GeV [[|l|. This lower limit also applies to the MSSM for small 
tan/?, even if squark mixing is maximal. In the CMSSM, maximal mixing is not attained, 
and the e + e~ — > Z° + h production rate is very similar to that in the Standard Model p2 



for all values of tan/3. As is well known, a 2.9-a signal for a Higgs boson weighing about 
115±o;| GeV has been reported gT[. 

In addition, the BNL E821 experiment has recently reported jy| a new value for the 
anomalous magnetic moment of the muon: g^ — 2 = 2 x a M , which yields an apparent 
discrepancy with the Standard Model prediction at the level of 2.6 a: 5a^ = (43 ± 16) x 
10~ 10 . This result has been well studied in the CMSSM [[0J, resulting in a strong preference 
to // > and relatively low values of m.1/2 and m . 

In addition to the phenomenological constraints, one must also carefully consider the 
resultant relic density of the LSP jOJ. This constraint generally disfavors large values of 
either of the SUSY breaking mass scales mo,mi/2 except in some well defined regions of 
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parameter space where an enhanced annihilation cross section ensures that the relic density 
flh 2 < 0.3 [46 ,4T| . These include regions where co-annihilations are important [48|, s-channel 



pseudo-scalar exchange is important [16], or in the focus point region at large m [Pfl . The 
co-annihilation region is important at large m 1( / 2 and the allowed regions due to pseudo-scalar 
exchange occur at large rai/2 and large mo at large tan/?. As we will see, the sensitivity to 
the phase 9 a occurs at relatively low mi/2 and mo. This is the region favored by the recent 
results from the g — 2 experiment. 

In Figure |], we show the cosmological and phenomeno logical constraints in the mi/2, m^ 
parameter plane for tan/5 = 10. The very dark (red) regions correspond to either m^ < m x , 
or mi x < m x where the subscript 1 denotes the lighter of the f or t mass eigenstates. These 
regions are ruled out by the requirement that the LSP be neutral. For A < 0, the excluded 
region is found at low values of mo where fi is the LSP. For Aq > 0, the fi LSP region 
is similar, but there is in addition a region where the stop is lighter than the neutralino 
at low mi/2 and extends to m ~ 500 GeV (in much of this region the stop is actually 
tachyonic). We show as (red) dash-dotted, nearly vertical lines the m^ = 113 GeV contour 
calculated using FeynHiggs f5Tfl . For A < 0, the limit is quite strong and excludes values 
of mi/2 ^ 440 GeV. For Aq > 0, the Higgs mass limit is not significant. The (dashed) 
bound on the chargino mass from LEP excludes very low values of m 1( / 2 < 140 GeV. The 
branching ratio for b — > sj excludes a dark (green) shaded area at low m 1( / 2 . For these 
parameter choices, this region is only present for A > 0. The (pink) medium-shaded region 
corresponds to the region preferredby the recent g—2 experiment. The region bounded by the 
solid curves corresponds to the 2-cr preferred region, whereas the dashed curves correspond 
to l-o". Finally, the (turquoise) light shaded region corresponds to the parameter space for 
which the neutralino relic density falls in the range 0.1 < Qh 2 < 0.3. 

The CP violating phases in the CMSSM have been studied 0-§J in relation to their 
effects on EDMs. Generally, it is found that while the phase of \x is strongly constrained by 
the EDMS to be close to either of its two CP conserving values, the phase of A is largely 
unconstrained. For this reason, in what follows, we will set 9^ = and concentrate on the 
effect of 9a on the rate for b — *■ sj. 

In the CMSSM, we specify the phases at the GUT scale. Since 9^, is not affected by 
the running of the renormalization group equations, the value of 9^ at the weak scale is 
identical to its input value. However, the real and imaginary parts of A run differently and 
hence, the weak scale phase, 9a differs from its input value 9a - In Figure ||, we show the 
resulting phase of A t at the weak scale as a function of the input value for several choices 
of m 1/2 = 200,400, and 600 GeV, with m = 400 GeV and |A | = 1000 GeV. The shaded 
regions show the dependence with respect to tan/3 between 10 and 30. The values of 9a 
for tan/5 = 10 are given by the lower edge of each shaded region, whereas the result for 
tan/5 = 30 is given by the upper edge. As one can see in the figure, the tendency of the 
RGE's is to drive A t towards real values, particularly at higher values of mi/ 2 . Nevertheless, 
the resulting phase at the weak scale is sufficient to have a strong impact on the branching 
ratio of b — > sj as we now discuss. 

In Figures |] - §, we show the calculated branching ratio of B — > A s 7 (scaled to the 
experimental value of 3.11 x 10 -4 ) as a function of m\/2 for 9a = 0, 7r/2, and 7r for both 
ji < (dotted curves) and // > (dashed curves). The horizontal lines correspond to the 95 
% CL range (2.27-4.01) xl0~ 4 . We have fixed the value of m = 140, for tan/5 = 10, which 
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the branching ratios for 9a = 0, n/2, and ir. Curves for /i > are shown as dotted and red, 
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experimental range. 

is the value which best yields a good relic density (for this value of \Aq\). As one can see the 
curves calculated for // < show a branching ratio which is too large unless m.1/2 ^ 430 GeV. 
At low mi/2, even though the phase of Aq has a strong effect on the calculated branching 
ratio, it can not reduce it to the experimentally allowed value. On the other hand, for /i > 0, 
we see that while all values of m.1/2 are allowed by b — > 57 for 9a = ft, there is a strong 
dependence on 9 a thus enabling one to set a limit on the phase of A from b — > 57, at least 
for some values of m.1/2 ■ One should bear in mind that while the range m.1/2 < 440 GeV is 
nominally excluded by the Higgs mass limit, the theoretical uncertainty in the calculated 
Higgs mass is about 3 GeV. The rrih = 110 GeV contour would lie at m.1/2 = 320 GeV in 
Figure £|b. 

In Figures and || we show the analogous behavior of the branching ratio for tan /3 = 20 
and 30 respectively. For the larger values of tan/3, it is not possible to satisfy the cosmological 
constraint with a single value of mo- However, we have chosen a value of mo which best fits 
the cosmological region for both signs of \x and A. In addition, we note that the dependence 
of the branching ratio on mo is relatively weak. Thus had we in fact varied mo with m.1/2 (to 
insure a proper relic density), the curves in Figures and [8] would differ only very slightly. 

The sensitivity of the branching ratio to the phase of A is also strong at the larger values 
of tan (3 as seen in Figures [7] and ||. As tan (3 is increased, one is pushed to higher values of 
mi/2 and therefore we are able to set bounds on 9a for a wider range in m.1/2. The Higgs 
mass bounds are also weaker at higher tan/3. 

As noted above, for a given range in m.1/2, we are able to use the constraints from b —> sj 
to set a limit on 9a - These limits are summarized in Figure |]. Shown there is the lower 



16 



3.5 



§" 2.5 

pa 

J 2 

pi 

5 1.5 



0.5 



e A = 7c/2\ \ e A = o 

= 7t 


m = 220GeV . 
IAqI = 1000 GeV 
tan P = 20 


--•'" 




- =71/2 .---"" ..--'" 

, .- 





100 150 200 250 300 350 400 450 500 

m 1/2 (GeV) 

FIG. 7. As in Fig. g for tan/? = 20 and m = 220 



3.5 



5F 2 - 5 

Pi 

Pi 
PQ 1.5 



0.5 







= 71/2 \ 




\ e,= 


= 






= 7C""- 












m = 

IA„I = 


320 GeV 
1000 GeV 


- 












tanP 


= 30 


= 71 


1 

.L - - - 


_ m 





i 


■■;::: 


'■"I-"-"* 


ffff: ::" = : = 




t 


^*'~ 


.--* 


,.-•' 








= 71/2 

■ * . - * 


>' 


..'' 


\- 


= 






" 



100 150 200 250 300 350 400 450 500 

m 1/2 (GeV) 

FIG. 8. As in Fig. g for tan/? = 30 and m = 320 



17 




150 



200 



250 



300 



350 



m 1/2 (GeV) 

FIG. 9. The lower bound on 9a as a function of mi/ 2 for tan/3 = 10,20 and 30. Values of 
Trio are the same as in Figures [g - [^. Note that at mj/2 ~ 170 GeV there will be a small window 
at 6a ~ which is allowed for tan/3 = 30. We do not include this curve in the figure. 

limit to 9a as a function of mi/2 f° r the three values of tan (3 indicated. Values of m were 
taken from Figures ^| - 18|. One can also read off from this figure the minimal value of vfi\i 2 
such that all values of the phase, 9a , are allowed. 

Finally, in Figures [10] and [11], we show the branching ratio as a function of tan/3 for fixed 
m = 400 GeV , \A \ = 1000 GeV, and m 1/2 = 200 and 400 GeV respectively. Recall that 
the sensitivity of these curves to mo (at least for relatively low mo) is rather small. As one 
can see, at small values of mi/2, there is a strong dependence of the branching ratio on the 
phase. At higher values of mm, not only is the dependence weaker, but also as one can see, 
for a given value of tan/3, either all phases are allowed or forbidden by the experimental 
constraint on b — > 57. 



V. SUMMARY 



In this work we have performed a thorough study of the constraints on the SUSY soft 
phases from B — > A s 7 decay using existing experimental bounds. Our analysis was restricted 
to the MFV scheme in which flavor violation occurs only through the CKM matrix. We 
considered both unconstrained as well as constrained supersymmetric models. 

In particular, we considered an unconstrained supersymmetric model, for which low 
energy supersymmetry bears no imprint of the stringy boundary conditions at ultra high 
energies. In particular, we chose a particle spectrum such that the charged Higgs boson, 
charginos, and the lighter stop are relatively light (with masses of order the weak scale). In 
this case, we found that: (i) the sum of the phases of the \x parameter and the stop trilinear 
coupling must take values around the CP-conserving point n, with a width determined by 
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experimental uncertainties; (ii) the sum of the phases of the \x parameter and the gluino 
mass tends to zero within the present experimental precision; (Hi) the LO CP asymmetry 
swings from —8% to +8% where positive and negative values are equally possible; (iv) 
the inclusion of the tan /3-enhanced corrections widens the allowed range of tan/3 values, 
though it does not lead to a significant change in the size of the CP asymmetry except for 
the fact that (a) it can be approximately twice as large as the LO prediction, and (b) it tends 
to follow the sign of the gluino mass for most of the parameter space, for tan/3 ~ O(50). 
Consequently, there exist observable effects of the SUSY threshold corrections at large tan (3. 
If experiments measure a large CP asymmetry (compared to the SM expectation), this may 
be interpreted as having originating from the soft phases and as an indication for weak scale 
SUSY. 

After analyzing the implications of the CP violating phases in an unconstrained su- 
persymmetric model, we turned to a detailed discussion of the CMSSM with explicit CP 
violation. In our numerical estimates, we took the // parameter to be real as implied by 
earlier studies of the EDM constraints, and we explored the b — > 57 branching ratio in re- 
gions of the parameter space allowed by the cosmological as well as other collider constraints 
displayed in Figure |j. In particular, we discussed the impact of the phase of the universal 
trilinear coupling at the GUT scale on the branching ratio. As the associated figures suggest, 
the branching ratio is quite sensitive to 6a , at least for moderately small values of m\/2- As 
we have shown, for positive values of //, while 9a = 7T, generally produced branching ratios 
in agreement with the experimental bounds, deviations from this CP conserving point can 
lead to branching ratios which are not compatible with experiment. We have also seen that 
the phase of the universal trilinear coupling can not ameliorate the inconsistency at low 7711/2 
when fj, < 0. 
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APPENDIX A: MASSES AND MIXINGS OF SUSY PARTICLES 

Given that the intergenerational mixings are proportional to the corresponding quark 

masses, the weak eigenstate squarks are approximately the mass eigenstates for the first two 

generations. However, the third generation squarks as well as the gauginos and Higgsinos 

can mix strongly after electroweak symmetry breaking. Since the chargino mass matrix M x ± 

— 2 + — 2 + 

is not hermitian, it is convenient to consider Ml = M ± • M x ± and Mr = M x ± ■ M ± which 

are hermitian and can be diagonalized by 

a m: Ca = aag. « , mi) with c. = ( .™% ~ sm J;C" ) ■ ( e 7 X ) < A1 > 

2 

The squark mass matrices M~, are similarly diagonalized. The relative sign between A t ^ 
and \x is set by the corresponding off-diagonal element of the sfermion mass matrix: m\ 2 = 
—m t (A t — //cot/3) for stops, and m\ 2 = —mb(Ab — //tan/3) for sbottoms. In eq. (|A.1|) , the 
explicit expressions for the angle parameters are given by 
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77= -Arg[A f * -//cot/3] , tan 2^ 



2m f |A* - //cot/3| 



M? - M? - (1/6) cos 2/3 (5Mf - 8M 2 ,) 



7r = -Arg [A* b - //tan/3] , tan 2% 



2m b |v4£ - //tan/3| 



'' JW? - M| + (1/6) cos 2/3 (Mf - 4M, 



ir) 



7 L = -Arg [M 2 + // cot /3] , tan 2fl L = 
7 R = -Arg [M 2 + //* tan/3] , tan 26> R 



^ sin /? |M 2 + // cot /?| 
M| + |/x| 2 + 2M^cos2/3 
v/8M w cos/3 |M 2 + //* tan/3| 
" Ml- |//| 2 -2M2 / cos2/3 



(A.2) 



where < ^ Q < 7r/2 to ensure M 2 > M 2 2 . Finally, one can choose r^ L = pj L = 0, and 



Vr = Ar £ 
Pr = Arj 



(A.3) 



Cfl (M 2 c L + V2M W sin/W^) + s fl e" i7 « (V2M W cos/3c L + //s^e*' 

c R {-V2M W cos/3s L e-^ - //c L ) - s R e^ R (-M 2 s L e-^ L + ^ sin/?c L )] (A.4) 



with sl^ = sin 6^^ and cl,r = cos^/j. The mixing matrix C a guarantees that both 
chargino masses are real positive with M 2 > M 2 and M~ > M~ . 



Finally, the neutralinos are described by a 4 x 4 mass matrix 



M u 



(Mi Mzs^cos^ -M z s w sin (5\ 

M 2 -M^c^, cos /3 M z c w sin /3 

M^s^cos/3 —Mzc w cos /3 

V — M^s^, sin /? M z c w sin (3 — // 







/ 



which can be diagonalized numerically via 

C£m°C = diag. (M x? , ■ • • , M x o 



where M v o < ■ ■ • < M v o . 



(A.5) 



(A.6) 



APPENDIX B: LOOP FUNCTIONS 

The two-point function entering the expressions of ebb,ts,tb (Q) * s 

x , y 



W [x, y] 



(1 - x) (x - y) 



lnx + 



;i - y) (y- x) 



\ny 



(B.l) 



with H [1, 1] = -1/2, and H [1, 0] = -1. 

The loop functions entering the expressions of the Wilson coefficients fl9HTTl) are given 
by 

. rri x (7 - 5x - 8x 2 ) x 2 (3x-2) 1 
F 7 LL [x] = — ^ tt^ + -^ tt/ hix 



36(x - l) 3 



6(x-l) 4 



F, LK [x] 



5 — 7x x (3x — 2) 



+ 



6(x-l) 2 3(x-l 



lnx 
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7-i 7" T r 1 X [Z -\- OX X ) OX 
Ft \x\ = — ; "^ ; -r 111 X 

8 [ J 12(x- l) 3 6(x-l) 4 

7- r T /-? r n ~~ l - i 

^W-if^-2(^TF 1 -- < a2 > 

In order to make numerical estimates, it is useful to know the values of these functions 
evaluated at x — 1: 

F 7 LL [l]=-^, Fi L [l] = -±, F 7 LR [l] = ^, 

F 8 LR [1]=\, F 7 LL [1] = -L , F 8 LL [1] = -I . (B.3) 

The two-loop EDMs, on the other hand, depend on the two-loop function 

F[x]= / ^ ^-= In — (B.4) 

Jo x — ZZ X 

where z = 1 — z with F [1] ~ —2.34. 
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